Nuclear deformation effect on the binding energies in heavy ions 

Y. S. Kozhedub\ O. V. Andreev\ V. M. Shabaev^'^, 1. I. Tupitsyn\ 
C. Brandau^, C. Kozhuharov-^, G. Plunien"', and T. Stohlker^'^ 
^ Department of Physics, 
St. Petersburg State University, 
Oulianovskaya 1, Petrodvorets, 
St. Petersburg 198504, Russia 
^ Max-Planck-Institut fur Physik Komplexer Systeme, 
Nothnitzer Str. 38, 
D-01187 Dresden, Germany 
^ Gesellschaft fur Schwerionenforschung (GSI), 
Planckstrasse 1, D-64291 Darmstadt, Germany 
^ Institut fiir Theoretische Physik, 
TU Dresden, Mommsenstrasse 13, 
D-01062 Dresden, Germany 
^ Physikalisches Institut, 
Ruprecht-Karls- Universitat Heidelberg, 
Philosophenweg 12, 
D-69120 Heidelberg, Germany 

(Dated:) 



1 



Abstract 

Nuclear deformation effects on the binding energies in heavy ions are investigated. Approximate formu- 
las for the nuclear-size correction and the isotope shift for deformed nuclei are derived. Combined with 
direct numerical evaluations, these formulas are employed to reanalyse experimental data on the nuclear- 
charge-distribution parameters in ^^^U and to revise the nuclear-size corrections to the binding energies in 
H- and Li-Uke ^^^U. As a result, the theoretical uncertainties for the ground-state Lamb shift in 238^91+ and 
for the 2pii2 — 2s transition energy in 238^89+ significantly reduced. The isotope shift of the 2pj — 2s 
transition energies for ^^^Nd^^"*" and ^^^Nd^^"*" is also evaluated including nuclear size and nuclear recoil 
effects within a full QED treatment. 

PACS numbers: 
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I. INTRODUCTION 



As is known (see, e.g., Refs yi |2|]), the finite-nuclear-size correction to the atomic energy 
levels is sufficiently well determined by the root-mean-square (rms) radius of the nucleus. Fol- 
lowing Franosch and Soff [[30, the uncertainty due to this correction was usually estimated by 
adding quadratically two errors, one obtained by varying the rms radius within its error bar and 
the other obtained by changing the model of the nuclear charge distribution from the Fermi to the 
homogeneously-charged- sphere model. This rather conservative estimate was sufficient in so far as 
the total theoretical uncertainty was mainly determined by other contributions. The recent progress 



made in calculations of higher-order QED and electron-correlation corrections [|4,j5|, a 



and 



the current status of precision experiments with heavy few-electron ions 
require, however, a more accurate consideration of the nuclear shape and deformation effects. 
Corresponding improvements are given in the present paper. 

The finite-nuclear-size correction is studied both numerically and analytically. Approximate 
analytical formulas for this effect are derived in the general case of a deformed nucleus. Special 
attention is paid to evaluation of the nuclear-size correction to the binding energies of H- and Li- 



m. The 



like uranium, where the most accurate experimental data were recently reported jlBL 
study performed in the paper is employed to revise the value of the nuclear rms charge radius 
for ^^'^U and to recalculate the corresponding correction to the binding energies. As a result, the 
theoretical accuracy of the ground-state Lamb shift in 238^91+ ^j^^ ^^yq 2pi/2 — 2s transition 
energy in ^38^89+ significantly improved. 

The isotope shift of the 2pj — 2s transition energies for the isotopes A=142 and A=150 of Li- 
like '^Nd^^+ is evaluated as a function of the difference 5(r^) of the nuclear mean-square charge 
radius. The calculation includes the nuclear size correction to the one-electron Dirac binding 
energy as well as the corresponding effect on the electron-correlation, Breit interaction, and QED 
contributions. The mass shift including the nonrelativistic, relativistic, and QED recoil effects is 
also evaluated. Combined with an estimate of the nuclear polarization effect on the binding energy, 
these data can be used to extract the 5 (r^) value from the corresponding experiment [1150. 

The relativistic unit system {h = c = m = 1) and the Heaviside charge unit (a = e^/47r, e < 
0) are employed throughout the paper. 
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II. FORMULATION 



The Coulomb interaction between an atomic electron and the nucleus is given by 

e2 1 

V{fe,fi,. . . ,rz) = -—y^T^ — , (1) 

Att ^-^ Ire — rA 

1=1 

where is the electron position, rj is the position of i-th proton, and the summation runs over all 
protons of the nucleus. If we neglect nuclear polarization effects, we can restrict our consideration 
of the operator V to a model space, where the nuclear states may differ from each other only by 
the projection of the total angular momentum on the laboratory Z axis. 

In what follows, we assume that the nuclear Hamiltonian can be separated into rotational and 
intrinsic parts, the nucleus is axially symmetric and has reflection symmetry with respect to the 
plane which is perpendicular to the axial- symmetry axis. With this assumption, the nuclear wave 
function can be written as [1^ 1^ 

\IMK) = ^^^{Dj,,,i<l>, e, ^)x'Ar') + {-iy~'DL^^,,{<l>, Q, ^)x'-Kir')^ (2) 
forKy^O and 



|/M0) = J^LLLdU^, e, vl/)x^(r') = -^YjMie, $ )Xo(^') (3) 

for if = 0, where / is the total nuclear angular momentum, M and K are its projections on the 
laboratory and the nuclear body-fixed Z axis, respectively, A denotes the other intrinsic quantum 
numbers, and (—1)"' must be considered as an operator defined by its action on the wave functions 
for given intrinsic angular momenta Jl^ 17]. Here and below the prime indicates variables taken 



in the nuclear coordinate frame and r' denotes the whole set of the internal nuclear coordinates. 
The Euler angles $, 6, \1' in the Wigner D-functions give the orientation of the intrinsic body-fixed 
system with respect to the laboratory frame. 

For a given internal nuclear state with K = 0, we should average the interaction operator 
V{fe,fi, . . . ,rz) with the internal nuclear wave function Xoi^')- obtain 

where the nuclear charge distribution p(r') is defined by 

z 

W) = (Xo'lE'^(^"-^1)l^o) (5) 

i=l 
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and r denotes the position of the f' vector in the laboratory coordinate frame. With the assumptions 
considered above, the density p(F) can be expanded in terms of spherical harmonics as 

p(r-^) = po{r')Yoo{n') + P2{r')Y2o{n') + p,{r')Y,o{fi') + ■■■ (6) 

with the multipole components 

pi{r) = / dnp{f)Yio{n) , (7) 



where n = r/r. Making use of the usual spherical harmonic expansion of |re — r| \ expression 
([4]) can be written as 



;— n »n— V e / 



/=0 m=~l 



X p{f')^^YiUn)Y:m{ne). (8) 

To integrate over the nuclear angular variables, we transform Yimin) in Eq. ([8]) to the nuclear 
coordinate frame 

I 



m'=-l 



We have 



{Xo\V\Xo) = J2^2k{re), (10) 



where 



k=0 



oo 

= / drr'pi{r)(-j^e{r-re) + -j^eir,-r] 



47r 





I 



- yUne)D'omi^,Q,^l!). (11) 



2/ + _ 

m=—i 



In the following, we restrict our calculations of nuclear size effects on atomic binding energies 
to even-A nuclei with total spin / = in the ground state. An extention to non-zero nuclear 
angular momenta (/ 7^ 0) can be performed in a similar way. In the case / = 0, the interaction 
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potential © must be averaged with the nuclear state \IMK) = |000): 

oo / 

/ rlrr^ nA r\ I ^—f^(r — r \ -I- - 



2 ^ oo t „ / ' ' 

v{f,) = (000|l^|000) = - ^ 5^ X] / drr''pi{r) i -^^Q{r - r,) + i^e(re - r] 



1=0 m=—l 'q 

X — Yim(ne) c/9sin0 / c/\E' 



' e 



21 + 1 







-^D'ooi^, e, mU^, 0, v[/)^do„(<i., e, VI/) 

vStt^ vStt^ 



oo 

o2 



/Att 

Here po is defined by Eq. ([7]) 



/drr'^poir) ( -6(r — re) H Of^e — t) ) . (12) 
\r re J 





po(r) = / dnp{r)Yooin) = — = / dnp{r). (13) 



In terms of the usual spherically- symmetric nuclear charge density 

p{r) = ^ j dnpir) (14) 



we obtain 



v(re) = —AnaZ [ drr'^p(r)( -Q(r — Ve) -\ Bfrg — r)), (15) 

J \r J 



where a is the fine structure constant. Thus, if we restrict our consideration to the case / = 0, 
the summation over / disappears and the interaction potential becomes spherically- symmetric. To 
calculate the energy shift due to the finite-nuclear- size effect one has to solve the Dirac equation 
with the potential v{r) given by Eq. (flSl) . 

For deformed nuclei the nuclear charge density is usually described by a modified Fermi model 

Pir) = — 7- -—, (16) 

1 + exp [{r — c)/a\ 

with P parameterization of nuclear deformation 

oo I 

C = Co(l + f^l-^Ylrn{e,f)) (17) 

1=1 m=—l 

consistent with the normalization condition: 



d^rp{r) = 1. 
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Assuming axial symmetry and considering only quadrupole and hexadecapole nuclear deforma- 
tion, the expression (fTTI) reduces to 

C = Co(l+/320>20 + /940no). (18) 

Before turning to the numerical evaluations of the nuclear-size correction for some ions of 
experimental interest, we also present approximate analytical formulas that explicitely take the 
nuclear deformation into account. 



ni. APPROXIMATE ANALYTICAL FORMULAS FOR THE ENERGY SHIFT 

According to the method of Ref. \^ the calculation of the one-electron finite-nucleus-size cor- 
rection AE for an arbitrary nuclear model can be reduced to the calculation of AE for the model 
of a homogeneously charged sphere with an effective radius R. To a high degree of accuracy, the 
effective nuclear radius for j = 1/2 states is given by M\ 



R 



1/2 



where 



(19) 



1 



drp{r)r^ = I drr^~^'^ p{r) . 



Then AE can be evaluated using the following approximative formulas [|20 



AE, 



ns 1 
7 



^— ^ 1+ «Z Vn.i 2— - ^ 
lun ^ \ nXc 



27 



mc 



{aZYn" 



AE„ 



1 



7 



40 



l + {aZffnpAo^Z))(2 



aZR\ 



27 



nXc J 



mc 



f{aZ) =bo + hi{aZ) + h2{aZf + h{aZf . 



(20) 



(21) 

(22) 
(23) 



Here n is the principal quantum number, Ac = h/mc, and 7 = yl — {aZy. The coefficients 
60 — ^3 for a number of states are given in Ref. [2]. Formulas (fT9l)-(l23l) allows one to calculate AE 
in the range Z = 1 — 100 with a relative accuracy of 0.2%. For more precise formulas we refer to 



Refs. U, 



1 



For the deformed Fermi distribution given by Eqs. (fT6l) . (fTSl) within the /3|o and /3|o approxi- 
'lo ~ /5|o) 

3 f 'Ka\^ 3/^ 3 ^A^o 1 /5 



mation (as a rule, /5|q ~ /5|q), we derive 



N 



(^20 + 4^*^40 



(24) 



5 

+ 2^ 



9 ^ / vra V A 3 „ 

1 + -^/?40 + — 1 + 77^/^40 



/^20 



+ 



5 



427r V TT 



/5: 



20 



112772^20 +27r 



7ra\ 
Co/ 



'40 T' 



(25) 



+ 



21 

An 



5 /5 
47r V vr 



+ 



75 
16^ 



l + 2(™) 

Co / 15 V Co 

2 



/540 
4n 



vra 
Co 



+ 



' A.) f - 

TT / V Co 



/^20 



20 



^20 + 4^ 



^ ^ 10 / 7ra\ _|_ 7 / 7ra\ 
3 V Co / 3 V Co / . 



40 



(26) 



Expanding (r^) and (r^) in terms of the /5 parameters and keeping the two lowest-order terms 
yields 



/^2\ _lfo 2 , 7„2 2N , 7cg + 3(7ra)2 3 



/52^o + ° 



9cn + (vra) 



■ CO ' ' ^ CO 

ir') =i(3c^ + ISvr^a^c^ + 31vrV) + + '^"'f j + ^/3,^o 

^ CO ' 



/^20' 



(27) 



+ 



17c^ + 32vrVc2 + 3vrV 3 



1 + (^)2 

V CO / 



98vr V vr 



_^20- 



(28) 



In the limit, where P20 tends to zero, the ordinary Fermi distribution is recovered. Substituting 
Eqs. (I27l)-(l28l) into formulas (fT9l)- (l22l) . one immediately finds AE for a hydrogenlike atom with a 
deformed nucleus, provided the parameters cq, a, and P20 are known. 

To study the role of nuclear deformation in calculations of the finite-nuclear-size correction, let 
us consider the energy difference for two isotopes. Since this difference can be approximated as 

m 



6E = AE2 - AEi ~ 2-f{SR/R)AE, 



(29) 
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we have to find the dependence of SR/ R on variations of the nuclear charge distribution parame- 
ters. Assuming that the value 1 /N, which determines the nuclear volume at a = 0, is proportional 
to atomic number A, we derive 



6R 16A l 7i^6{a^) ^.,^2 



(30) 



where the first term is due to an increase of the nuclear volume, the second one results from 

If the 



a change of the parameter a, and the third one represents nuclear deformation 11191 . 
spherically-symmetric nucleus is considered as a reference (5(/9|o) = /^lo)' ^^'^ the parameter a is 
the same for both isotopes, we get 

6R 16A 5 ^2 

— - — 7- + —P2()- (31) 

This formula gives a simple way to determine the nuclear deformation parameter P20, provided the 
isotope shift is known, e.g., from experiment. 

Alternatively, considering (r^)^/^, a, and /32o as free independent parameters, we obtain 

3 



5R 5irY'^ 3 , ^.o^r'^fa') 

— ~ — ^ — (aZY — - 

R ~ (r^y/^ 70^ ' (r2) 



567r 



(32) 



(r2)i/2 

This formula shows that, to a good accuracy, the isotope shift is determined by the change of the 
rms radius. 



IV. NUCLEAR SIZE CORRECTION TO THE BINDING ENERGIES IN 238u9i+ ^ND 238u89+ 



In this section the formulation given above is applied to deduce a new value for the rms radius 
in ^^^U and, with this value, to revise theoretical predictions for the ground state Lamb shift in 
238^91+ j-j^g 2pi/2 — 2s transition energy in 238^89+ 



Compilation of the rms values [|21 



22 



2311 employed experimental data for nuclear charge 
distribution parameters obtained by various experimental methods. In case of 235, 238^ most 
recent compilation by Angeli [2l[ 22] includes data from elastic electron scattering [24], muonic 



atom X-rays [|25L |26D, X-ray isotope shifts 



ifts 



2711 . and optical isotope shifts pSQ . Since in 



Ref. [25] the experimental data are given in terms of the parameters a, cq, (^20, and P^o, one should 



2211 this was achieved based on formulas 



first evaluate the corresponding rms values. In Refs. 1|21L 
which only partly account for the deformation effect. In the present work we improved the Angeli's 
evaluation employing formulas (|24l) - (|25l) as well as the direct numerical calculations. As a result. 
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260 



we obtained the (r^)^/^ values which are close enough to those from the other sources 1124 . 
In case of ^^^U, the compillation of the improved data for (r^)^/^ and the 5(r^) data from Refs. 
H 0,13 28], performed by Angeli [29], yields (rY^^ = 5.8569(33) fm. This value differs 
from the corresponding value from the previous compillation, (r^)^/^ = 5.8507(72) fm [21]. As 
to the other nuclear-charge-distribution parameters, in accordance with the available experimental 
data 125, Q, we use a = 0.50(5) fm, /^so = 0.27(1), and (3^0 = 0.05(10) assuming rather 
conservative errors bars. These parameters differ from those employed in similar calculations by 



Blundell et al. l|30 



X-ray experiment 112611 



and by Ynnerman et al 



. iH, 



who adopted exclusively the data of the muonic 



The finite-nuclear- size correction is obtained by solving the Dirac equation with the potential 
([T5l) and taking the difference between the energies for the extended and the point-charge nu- 
cleus. In order to investigate the importance of the nuclear deformation effect, the calculations of 
the finite-nuclear-size correction are also performed using a spherically-symmetric nuclear charge 
distribution with the same rms value or with the same nuclear volume. The results of these calcula- 
tions are compared with each other in Table HJ In addition to the direct numerical (N) calculations, 
the analytical (A) results obtained by formulas ([T9l)- (|26l) . which provide a 0.2% accuracy, are pre- 
sented as well. As one can see from the table, if the rms value is kept to be the same, the nuclear 
deformation provides a 0.06% energy shift. If the nuclear volume is constant, the energy shift 
amounts to about 2%. It can also be seen that the energy shifts obtained by analytical formulas 
(|29l)-(l32]) are in a reasonable agreement with the exact numerical results. We note also that the 
effect of hexadecapole deformation (~ P^q) is extremely small for ^^^U, provided the rms radius 
is kept to be constant. 

Thus to calculate the nuclear size correction for 238y9i+ ^ 0.1% accuracy one needs to 
account for the nuclear deformation effect. Finally, the nuclear-size corrections for 238^91+ 



are AE(ls) = 198.54(19) eV, AE{2p 



1/2 



2s) 



-33.304(30) eV, and AE{2py2 - 2s) 



-37.714(34) eV. 

In the last compilations of the ground-state Lamb shift in 238^91+ ^ ^ ^^le 2pi/2 — 2s 
transition energy in 238^89+ ^ gjj [q^^x theoretical uncertainties were mainly determined by the 
finite-nuclear-size corrections. The new values for these corrections obtained in the present work 
provide significant improvements of the theoretical predictions for both H- and Li-like uranium. 
In Table HI] we present individual contributions to the Is Lamb shift in 238^91+ jj^g uncertainty 
of the total theoretical value, 463.99(39) eV, is now mainly determined by uncalculated two-loop 
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Table I: The exact numerical (N) and approximate analytical (A) results for the finite-nuclear-size correction to the energies of Is, 2s, and 2jOi/2 states of 
238u9i+ ((^2^1/2 ^ 5.8569(33) fm, a = 0.50(5) fm, /^ao = 0.27(1), and /?4o = 0.05(10)), in eV. The results for a deformed (D) nucleus are compared 
with the results obtained for a spherically-symmetric nuclear model with (1) the same value of the rms value ((r2)V2 = {r'^)]^'^) or with (2) the same 
nuclear volume (l/AT = {1/N)-d). 



Nuclear 


(r2)V2 




^20 


/340 


Co 


Is 


2s 


2pi/2 


2pi/2 


Method 


model 


















-2s 






(fm) 


(fm) (fm) 






(fm) 


(eV) 


(eV) 


(eV) 


(eV) 




Def. nuc. 


5.8569(33) 6.2384 0.50(5) 0.27(1) 0.05(10) 7.0140 198.54(19) 37.714(34) 4.410(4) 


-33.304(30) 


N 














198.39 


37.651 


4.412 


-33.239 


A 


(1) Sph. sym. 


5.8569 


6.2088 0.50 


0.00 


0.00 


7.1704 


198.68 


37.740 


4.413 


-33.327 


N 














198.61 


37.692 


4.417 


-33.275 


A 


(2) Sph. sym. 


5.7805 


6.1303 0.50 


0.00 


0.00 


7.0663 


194.90 


37.025 


4.328 


-32.696 


N 


l/N = (l/iV)D 












194.77 


36.963 


4.331 


-32.632 


A 



Table II: Individual contributions to the ground-state Lamb shift in 238^91+^ 



Contribution 
Finite nuclear size 
First-order QED 
Second-order QED 
Nuclear recoil 



Value Reference 
198.54(19) This work 
266.45 [321 
-1.26(33) [4J 
0.46 [3J] 



Nuclear polarization -0.20( 10) [34, 35J 



Total theory 
Experiment 



463.99(39) 
460.2(4.6) 



[13] 



QED corrections, in particular, the mixed vacuum-polarization self-energy contribution jSq] . The 
obtained result is in a good agreement with the recent experiment lll3ll . 

TableUnipresents individual contributions to the 2pi/2 — 2s transition energy in 238^89+ Qq^_ 
pared to Refs. M,M, it contains the new value for the nuclear-size correction and the new value for 
the three- and more photon effects. The latter correction was evaluated within the Breit approx- 
imation employing the large-scale configuration-interaction Dirac-Fock-Sturm (CI-DFS) method 



1411 



421] . The procedure successfully used for Li-like scandium 14311 was applied here as wel 



. For 



uranium, we report a good agreement with the previous evaluations of this correction 144 , 
The uncertainty ascribed to this correction incorporates all three- and more photon effects which 
are beyond the Breit approximation. The entry labeled "Screened QED" represents the sum of the 
lowest-order self-energy and vacuum-polarization screening diagrams Issi 3^ . 

Table Uni shows that now, after our revision of the finite-nuclear-size correction, the total theo- 
retical uncertainty is mainly influenced by higher-order QED effects. The total theoretical value 
of the transition energy, 280.71(10) eV, agrees well with the most precise experimental value, 
280.645(15) eV lll4ll . Comparing the first- and second-order QED contributions with the total the- 
oretical uncertainty, we conclude that the present status of the theory and experiment for Li-like 
uranium provides a test of QED on a 0.2% level to first order in a and on a 6.5% level to second 
order in a. 
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Table III: Individual contributions to the l-piji — 2s transition energy in ^■^^'u , in eV. 



Contribution 


Value 


Reference 


One-electron nuclear size 


-DD.D\J\D ) 


This work 


One-photon exchange 




This work 


One-electron first-order QED 




1 j2I 


Two-photon exchange within the Breit approx. 


-13. j4 


\il 1 


Two-photon exchange beyond the Breit approx. 


U. 1 / 


1 J / 1 


Screened QED 


l.lu 




One-electron second-order QED 


V.ZZ{p) 




i nree- aiiQ inore pnoion enccis 


V. l^y 1 ) 


inis worK 








Nuclear polarization 


0.03(1) 


[34,35J 


Total theory 


280.71(10) 




Experiment 


280.645(15) 




Experiment 


280.59(10) 


m 


Experiment 


280.52(10) 


[12] 



V. ISOTOPE SHIFT OF THE 2pj - 2s TRANSITION ENERGIES FOR 42^^07+ ^ND ^^^Nd 



i57+ 



150iv^57+ 



In this section we evaluate the isotope shift of the 2pj — 2s transition energies for the iso- 
topes A= 142 and A= 150 of Li-like ^Nd^^+, where the '^"Nd nucleus is strongly deformed 
{(320 = 0.28(5), see, e.g., Ref. MM)- To date, there are about 20 publications, where the mean 



21 



square charge radius difference 5(r^) for these isotopes is reported (see Refs. Iil5 . 
and references therein). Apart from some outliers, the majority of the experimental data cover a 
range from about 150,142^^^2^ ^ j^g 150,142^^^2^ ^ gg ^^^^2 p^j. j-j^j^ reason, we perform 

calculations of the isotope shift for the entire range of ^^°'^'^^5(r^), from 1.18 to 1.38 fm^. With 
these data, one can easily find the value of 150,142^^^2^ from the experimental value of the isotope 



shift 



The isotope shift is given by the sum of the field shift, which is due to the finite-nuclear-size 
effect, and the mass shift, which is determined by the nuclear recoil effect. To evaluate the field 
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shift we used the large-scale CI-DFS method 114 iL l42h . with the Breit interaction included. The 
spherically-symmetric ^"^^Nd nucleus served as a reference with the rms 

radius of (r2) 1/2 = 4.9118 

fm from the compillation by Angeli [12 ill . The other nuclear parameters are taken to be a = 0.52(2) 
fm for both isotopes, /^ao = for ^^^Nd and /^ao = 0.28(5) for ^^^Nd [23]. We note that variations 
of these parameters within their error bars do not affect the isotope shifts at the accuracy level 
considered. 

The full relativistic theory of the nuclear recoil effect can be formulated only in the framework 



of QED [47]. To eva 



uate the recoil effect within the lowest-order relativistic approximation one 



can use the operator 148, 

Hm '■ 



m 



-E 



2M 



aZ 



OLi 



[OLi ■ ri Vi 



■Pi 



(33) 



where M is the nuclear mass and Pi is the momentum operator acting on the i-th electron. The 
expectation value of Hm on the many-electron wave function of the system, obtained by the 
CI-DFS method, yields the recoil correction to the energy levels to all orders ml/ Z within the 
{aZYm? / M approximation. The recoil correction which is beyond this approximation is termed 
as the QED recoil effect. For the 2pi/2 — 2s (2^3/2 — 2s) transition the mass shift comprises of 
2.44 meV (2.53 meV) from averaging the nonrelativistic part of the recoil operator (the first term 
in Eq. ((331) ) with the relativistic many-electron wave function, —1.14 meV (—1.03 meV) from the 
relativistic par t (the second term in Eq. ([SSl)). and of 0.33 meV (0.30 meV) from the QED recoil 
effect [|4o[ Isoll . The recoil correction of the next order in m/M is negligible in the case under 
consideration. Finally, the total mass shift sums up to 1.63 meV for the 2]9i/2 — 2s transition and 
to 1.80 meV for the 2^3/2 — 2s transition. 

Next, one should account for the influence of the nuclear size variation on the one-loop QED 
corrections. Using comprehensive tabulations for the nuclear-size correction to the self-energy 
contribution [51] and evaluating the corresponding effect on the vacuum-polarization contribution, 
we derive 0.2 meV for the QED correction to the isotope shifts under consideration. 

Finally, we have to consider the nuclear polarization (NP) effect. This correction is determined 
by the electron-nucleus interaction diagrams in which the intermediate states of the nucleus are 



3511 ■ Since the NP correction 



excited. This effect was evaluated for a number of ions in Refs. [|34l . 
is most sizeable for deformed nuclei, we estimated it for ^^^Nd taking into account the transition 
to the first excited 2+ state at 130.21 keV only. Taking the nuclear transition probability from Ref. 
ll52ll and evaluating the sum over intermediate electron states numerically as well as analytically 
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Table IV: Isotope shift for the 2pij2 — 2s transition in Li-Uke 150,142^^^57+^ rpj^g ^^^^ ^j^j^^ includes 

one-electron Dirac, electron-correlation, and Breit-interaction contributions. The mass shift incorporates 
nonrelativistic, relativistic, and QED recoil effects. The QED correction represents the sum of one-loop 
self-energy and vacuum-polarization contributions. 



6{r^) 


Field shift Mass shift 


QED 


Nuc. pol. 


Total 


(fm2) 












1.180 


-0.0366 


0.0016 


0.0002 


0.0003 


-0.0345 


1.200 


-0.0372 


0.0016 


0.0002 


0.0003 


-0.0351 


1.220 


-0.0379 


0.0016 


0.0002 


0.0003 


-0.0358 


1.240 


-0.0385 


0.0016 


0.0002 


0.0003 


-0.0364 


1.260 


-0.0391 


0.0016 


0.0002 


0.0003 


-0.0370 


1.280 


-0.0397 


0.0016 


0.0002 


0.0003 


-0.0376 


1.300 


-0.0403 


0.0016 


0.0002 


0.0003 


-0.0382 


1.320 


-0.0410 


0.0016 


0.0002 


0.0003 


-0.0389 


1.340 


-0.0416 


0.0016 


0.0002 


0.0003 


-0.0395 


1.360 


-0.0422 


0.0016 


0.0002 


0.0003 


-0.0401 


1.380 


-0.0428 


0.0016 


0.0002 


0.0003 


-0.0407 



according to formulas derived in Ref. pSl] . we obtain 0.3(3) meV for the nuclear polarization 
contribution to the isotope shift for both transitions. 

The results of our calculations are presented in TablesHVlandrVlfor the 2pi/2 — 2s and 2p3/2 — 2s 
transitions, respectively. With the numbers compiled in these tables, one can easily deduce the 
nuclear mean-square charge difference 5(r^), provided the isotope shift is known from experiment 
lIlSll . In addition, using formula (|3T1) one can derive the quadrupole deformation parameter /?2o to 
an accuracy of about 20 — 30%. 



VI. CONCLUSION 

The finite-nuclear-size correction to the binding energies in heavy ions has been studied in this 
paper. In the general case of a deformed nucleus, approximate analytical formulas for this effect 
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Table V: Isotope shift for the 2^3/2 — 2s transition in Li-like I5u,i4ij^j5/+^ jj^g ^^^^ ^^^^^ includes 

one-electron Dirac, electron-correlation, and Breit-interaction contributions. The mass shift incorporates 
nonrelativistic, relativistic, and QED recoil effects. The QED correction presents the sum of one-loop self- 
energy and vacuum-polarization contributions. 



6{r^) 


Field shift Mass shift 


QED 


Nuc. pol. 


Total 


(fm2) 












1.180 


-0.0379 


0.0018 


0.0002 


0.0003 


-0.0353 


1.200 


U.UJOJ 


U.UU i 


0.0002 


0.0003 


-0.0362 


1.220 


-0.0392 


0.0018 


0.0002 


0.0003 


-0.0369 


1.240 


-0.0398 


0.0018 


0.0002 


0.0003 


-0.0375 


1.260 


-0.0404 


0.0018 


0.0002 


0.0003 


-0.0381 


1.280 


-0.0411 


0.0018 


0.0002 


0.0003 


-0.0388 


1.300 


-0.0417 


0.0018 


0.0002 


0.0003 


-0.0394 


1.320 


-0.0424 


0.0018 


0.0002 


0.0003 


-0.0401 


1.340 


-0.0431 


0.0018 


0.0002 


0.0003 


-0.0408 


1.360 


-0.0437 


0.0018 


0.0002 


0.0003 


-0.0414 


1.380 


-0.0443 


0.0018 


0.0002 


0.0003 


-0.0420 



have been derived and direct numerical calculations have been performed. In the special case 
of ^^^U the study has been employed to revise the nuclear-charge-distribution parameters and to 
recalculate the binding energies in H- and Li-like uranium. As the result, the largest theoretical 
uncertainties for the ground-state Lamb shift in 238^91+ ^^yq 2pi/2 — 2s transition energy 

in 238^89+ j^^yg been removed. Now the total theoretical accuracy is mainly restricted by higher- 
order QED effects. Tables Ull and Hill demonstrate that our theoretical results agree well within the 
error bars with the most precise experimental data. 

We have also evaluated the isotope shift of the 2pj — 2s transition energies for i''2]sjjj57+ 
i50]vj^57+ ^Qj. (different values of the mean-square nuclear charge difference 5(r^). The calcula- 
tion of the field shift takes into account electron-correlation, Breit-interaction, and QED effects. 
The mass shift is evaluated within a full QED treatment. The nuclear-polarization correction is 
also estimated. The data obtained allow one to extract the 5(r^) value from the corresponding 
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experiment. 



Acknowledgments 

Valuable communications with I. Angeli are gratefully acknowledged. This work was sup- 
ported in part by INTAS (Grant No. 06-1000012-8881), RFBR (Grant No. 07-02-00126), DFG, 
DAAD, and GSI. Y.S.K. and O.V.A. acknowledge the support by the "Dynasty" foundation. 



[1] W. R. Johnson and G. Soff, At. Data Nucl. Data Tabl. 33, 405 (1985). 

[2] V. M. Shabaev, J. Phys. B 26, 1103 (1993). 

[3] T. Franosch and G. Soff, Z. Phys. D 18, 219 (1991). 

[4] V. A. Yerokhin, R IndeUcato, and V. M. Shabaev, Phys. Rev. Lett. 91, 073001 (2003). 

[5] A.N. Artemyev, V.M. Shabaev, V.A. Yerokhin, G. Plunien, and G. Soff, Phys. Rev. A 71, 062104 

(2005) . 

[6] V. M. Shabaev, O. V. Andreev, A. N. Artemyev, S. S. Baturin, A. A. Elizarov, Y. S. Kozhedub, N. 
S. Oreshkina, 1. 1. Tupitsyn, V. A. Yerokhin and O. M. Zherebtsov, Int. J. Mass Spectrom. 251, 109 

(2006) . 

[7] V. A. Yerokhin, P Indehcato, and V. M. Shabaev, Phys. Rev. Lett. 97, 253004 (2006). 

[8] V. M. Shabaev, A. N. Artemyev, D. A. Glazov, 1. 1. Tupitsyn, A. V. Volotka, and V. A. Yerokhin, AIP 

Conference Proceedings, 869, 52 (2007). 
[9] J. Schweppe, A. Belkacem, L. Blumenfeld, N. Claytor, B. Feinberg, H. Gould, V. E. Kostroun, L. 
Levy, S. Misawa, J. R. Mowat, and M. H. Prior, Phys. Rev. Lett. 66, 1434 (1991). 
[10] S. R. Elliott, R Beiersdorfer, and M. H. Chen, Phys. Rev. Lett. 76, 1031 (1996); 77, 4278 (1996). 
[11] S. R. Elliott, P. Beiersdorfer, M. H. Chen, V. Decaux, and D. A. Knapp, Phys. Rev. C 57, 583 (1998). 
[12] C. Brandau, C. Kozhuharov, A. MuUer, W. Shi, S. Schippers, T. Bartsch, S. Bohm, C. Bohme, A. 
Hoffknecht, H. Knopp, N. Griin, W. Scheid, T. Steih, F. Bosch, B. Franzke, R H. Mokler, F. Nolden, 
M. Steck, T. Stohlker, and Z. Stachura, Phys. Rev. Lett. 91, 073202 (2003). 
[13] A. Gumberidze, T. Stohlker, D. Bana, K. Beckert, P Beller, H. F. Beyer, F. Bosch, S. Hagmann, C. 
Kozhuharov, D. Liesen, F. Nolden, X. Ma, P. H. Mokler, M. Steck, D. Sierpowski, and S. Tashenov, 
Phys. Rev. Lett. 94, 223001 (2005). 



17 



[14] P. Beiersdorfer, H. Chen, D. B. Thorn, and E. Trabert, Phys. Rev. Lett. 95, 233003 (2005). 

[15] C. Brandau, C. Kozhuharov, Z. Harman, A. Miiller, S. Schippers, Y. S. Kozhedub, D. Bernhardt, S. 
Bohm, J. Jacobi, E. W. Schmidt, F. Bosch, H. J. Kluge, P H. Mokler, T. Stohlker, K. Beckert, P 
Beller, F. Nolden, M. Steck, A. Gumberidze, R. Reuschl, U. Spillmann, F. J. Currell, 1. 1. Tupitsyn, V. 
M. Shabaev, U. D. Jentschura, C. H. Keitel, A. Wolf, and Z. Stachura, to be published. 

[16] A. G. Sitenko and V. K. Tartakovskii, Lectures on the Theory of Nucleus (Moscow: Atomizdat, 1972). 

[17] P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer, New York, 1980). 

[18] V. M. Shabaev, Opt. Spectrosc. 56, 244 (1984). 

[19] L. Wilets, D. L. Hill, and K. W. Ford, Phys. Rev. 91, 1488 (1953). 

[20] C. S. Wu and L. Wilets, Annu. Rev. Nucl. Sci. 19, 527 (1969). 

[21] I. Angeh, At. Data Nucl. Data Tables 87, 185 (2004). 

[22] I. Angeli, Nuclear Data Services of the International Atomic Energy Agency. INDC(HUN)-033, 

September 1999, http://www-nds.iaea.org/indc_sel.html. 
[23] G. Fricke, C. Bernhardt, K. Heilig, L. A. Schaller, L. Schellenberg, E. B. Shera, and C. W. De Jager, 

At. Data Nucl. Data Tables 60, 177 (1995). 
[24] C. W. Creswell, Ph.D. thesis, Mass. Inst. Techn., unpubl.; H. deVries, C. W. deJager, and C. deVries, 

At. Data. Nucl. Data. Tables 36, 495 (1987). 
[25] D. A. Close, J. J. Malanfy, and J. P Davidson, Phys. Rev. C 17, 1433 (1978). 

[26] J. D. Zumbro, E. B. Shera, Y. Tanaka, C. E. Bemis, R. A. Naumann, M. V. Hoehn, W. Reuter, and R. 

M. Steffen, Phys. Rev. Lett. 53, 1888 (1984). 
[27] R. T. Brockmeier, F. Boehm, and E. N. Hatch, Phys. Rev. Lett. 15, 132 (1965). 
[28] A. Anastassov, Yu. P. Gangrsky, K. P. Marinova, B. N. Markov, B. K. Kuldjanov, and S. G. Zemlyanoi, 

Hyperfine Interact. 74, 31 (1992). 
[29] I. Angeli, private communication. 

[30] S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. Rev. A 41, 1698 (1990). 

[31] A. Ynnerman, J. James, I. Lindgren, H. Persson, and S. Salomonson, Phys. Rev. A 50, 4671 (1994). 

[32] P J. Mohr, G. Plunien, and G. Soff, Phys. Rep. 293, 227 (1998). 

[33] V. M. Shabaev A. N. Artemyev, T. Beier, G. Plunien, V. A. Yerokhin, and G. Soff, Phys. Rev. A 57, 
4235 (1998). 

[34] G. Plunien and G. Soff, Phys. Rev. A 51, 1119 (1995); 53, 4614 (1996). 

[35] A. V. Nefiodov, L. N. Labzowsky, G. Plunien, and G. Soff, Phys. Lett. A 222, 227 (1996). 

18 



[36] S. Zschocke, G. Plunien, and G. Soff, Eur. Phys. J. D 19, 147 (2002). 

[37] V. A. Yerokhin, A. N. Artemyev, V. M. Shabaev, M. M. Sysak, O. M. Zherebtsov, and G. Soff, Phys. 

Rev. Lett. 85, 4699 (2000); Phys. Rev. A 64, 032109 (2001). 
[38] A. N. Artemyev, T. Beier, G. Plunien, V. M. Shabaev, G. Soff, and V. A. Yerokhin, Phys. Rev A 60, 

45 (1999). 

[39] V. A. Yerokhin, A. N. Artemyev, T. Beier, G. Plunien, V. M. Shabaev, and G. Soff, Phys. Rev. A 60, 
3522 (1999). 

[40] A. N. Artemyev, V. M. Shabaev, and V. A. Yerokhin, Phys. Rev. A 52, 1884 (1995). 

[41] 1. 1. Tupitsyn, V. M. Shabaev, J. R. Crespo Lopez-Urrutia, 1. Draganic, R. Soria Orts, and J. Ullrich, 

Phys. Rev. A 68, 022511 (2003). 
[42] 1. 1. Tupitsyn, A. V. Volotka, D. A. Glazov, V. M. Shabaev, G. Plunien, J. R. Crespo Lopez-Urrutia, A. 

Lapierre, and J. Ullrich, Phys. Rev. A 72, 062503 (2005). 
[43] Y. S. Kozhedub, D. A. Glazov, A. N. Artemyev, N. S. Oreshkina, V. M. Shabaev, 1. 1. Tupitsyn, A. V. 

Volotka, and G. Plunien, Phys. Rev. A 76, 012511 (2007). 
[44] O. M. Zherebtsov, V. M. Shabaev, and V. A. Yerokhin, Phys. Lett. A 277, 227 (2000). 
[45] O. Yu. Andreev, L. N. Labzowsky, G. Plunien, and G. Soff, Phys. Rev. A 64, 042513 (2001). 
[46] V. A. Yerokhin, A. N. Artemyev, and V. M. Shabaev, Phys. Rev. A 75, 062501 (2007). 
[47] V. M. Shabaev, Phys. Rev. A 57, 59 (1998); Phys. Rep. 356, 1 19 (2002). 

[48] V. M. Shabaev, Teor. Mat. Fiz. 63, 394 (1985) [Theor. Math. Phys. 63, 588 (1985)]; Yad. Fiz. 47, 107 

(1988) [Sov. J. Nucl. Phys. 47, 69 (1988)]. 
[49] C. W. Palmer, J. Phys. B 20, 5987 (1987). 

[50] A. N. Artemyev, V. M. Shabaev, and V. A. Yerokhin, J. Phys. B 28, 5201 (1995). 

[51] T. Beier, R J. Mohr, H. Persson, and G. Soff, Phys. Rev. A 58, 954 (1998). 

[52] S. Raman, C. W. Nestor, Jr., and P Tikkanen, At. Data Nucl. Data Tabl. 78, 1 (2001). 



19 



